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L». I Abstract: First we present a theoretical framework to compute the fully differential cross 

^^ ' sections for the top-quark productions and their subsequent decays at hadron colliders, in- 

t~^ ■ corporating the bound-state effects which are important in the ti threshold region. We 

include the bound-state effects such that the cross sections are correct in the LO approxi- 
mation both in the threshold and high-energy regions. Then, based on this frainework we 
compute various kinematical distributions of top quarks as well as of their decay products 
at the LHC, by means of Monte-Carlo event-generation. These are compared with the cor- 



H \ responding predictions based on conventional perturbative QCD. In particular, we find a 

characteristic bound-state effect on the {hW~^)-{hW~) double-invariant-mass distribution, 
which is deformed to the lower invariant-mass side in a correlated manner. 
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1. Introduction 

At the CERN Large Hadron Collider (LHC), the top quark will be produced copiously. The 
cross section for the top-quark pair-production amounts to several hundred pb [j^, |2|, ^, and 
order 10^ top-quark events will be observed each year if the LHC runs with 14 TeV collision 
energy and achieves the designed luminosity. Collecting these top-quark events, detailed 
analysis on the properties of the top quark will be possible, such as precise determinations 
of its mass and width, structure of electroweak and strong interactions, and its spin proper- 
ties . The current world average of the top-quark mass measurements from the combined 
analysis of CDF and DO collaborations at the Fermilab Tevatron reads mt = 173.1 it 1.3 
GeV [^] (see also p). Furthermore, the top-quark production process is considered as a 
standard candle process. Namely, it serves understanding detector performances, e.g. jet 
energy calibrations, from comparisons of experimental measurements with theoretically re- 
liable or well-controllable predictions, for observables including jet topologies, backgrounds 
and underlying events. 

There have been many studies on top-quark production processes at the LHC* Update 
analyses on the total pair-production cross-section are presented in fTl H B , including the 



*See e.g. M, hq] for more complete review. 



next-to-leading-order (NLO) correction 0, |8[ and resummation of threshold logarithms ||9|, 
IC, ^ in QCD. Differential distributions including decays of the top-quark have also been 



known up to NLO |12, |T^, |^], and various distributions are investigated in |15, 1£, 17, 18 

Recently, tt invariant-mass (m^f) distribution near threshold has been investigated 
incorporating the bound-state effects [^, ^. The effects are found to be significant at 
the LHC, since (in contrast to the Tevatron) the gluon- fusion channel dominates the cross 
section and there are significant contributions from (the remnant of) the color-singlet tt 
resonances. 

In this paper we compute the fully differential cross sections for the top-quark pair 
productions and their subsequent decays at the LHC. In particular, we incorporate the 
bound-state effects, which are important in the tt threshold region, into the cross sec- 



tions. We extend the studies of [^, 21 1 and present a theoretical frame to incorporate the 
bound-state effects to the differential cross sections. Using this result, we compute various 
kinematical distributions of the top quarks and their decay products at hadron colliders, 
by developing a Monte-Carlo (MC) event-generator incorporating the bound-state effects. 
(There exist similar MC event-generators for computing the top-quark cross-sections in the 
tt threshold region at future e^e~ colliders [^ g3|, |4|.) Through the analysis, we eluci- 
date the nature of the bound-state effects at various stages: at the partonic matrix-element 
level, both with and without including the decay of the top quark, and in the kinematical 
distributions after incorporating the initial-state radiation (ISR) effects. Theoretically the 
fully differential cross sections contain more information on the bound-state effects than 
just the tt invariant-mass distribution; for instance, it is known that the top momentum 
distribution is sensitive to the resonance wave functions in momentum space |25| , |2q| . From 
a practical point of view, the differential cross sections are useful for studying effects of 
various kinematical cuts, detector acceptance corrections, detector calibrations, etc. 

The method for incorporating tt bound-state effects has been developed mainly in the 
studies of ti productions in e~^e~ collisions p7| , p^ , p9|. Formally, in the limit where we 
neglect the top-quark width, Ft — )• 0, bound-state effects can be incorporated by resumma- 
tion of the Coulomb singularities (asZ/S)", where f3 is the velocity of the top quark in the 
tt cm. frame. In contrast to the e~^e~ collision, at hadron colliders, ti pairs are produced 
in both color-singlet and octet states, and the (partonic) collision energy is not fixed. Due 
to the latter reason, we have to set up a theoretical framework which is valid both in the 
threshold region (m-tj — 2?7ij) and in the high-energy region (m-tj ^ 2?TT,t). The former re- 
gion is where the bound-state effects (Coulomb corrections) become significant and where 
the non-relativistic approximation is valid. On the other hand, in the latter region, the 
bound-state effects are not significant and the top quarks are relativistic. We present a 
framework which takes into account all the leading-order (LO) corrections in both regions. 
Namely, we incorporate all the {as//3)"' terms in the threshold region, while we include 
all the /S" terms in the relativistic region. (Some of the important subleading corrections 
are also incorporated.) Furthermore, we interpolate the two regions smoothly in a natural 
way. 

Another important aspect in computing the differential cross sections for the top-quark 



productions and decays is to construct full amplitudes corresponding to the bW^bW^ final 
state and to incorporate off-shellness of the top quarks appropriately ^^ . The former is 
important to incorporate the polarizations of t and i and angular correlations in their decay 
products. (In Appendix D, we will further discuss the polarizations of the M^'s from the 
top quarks.) The latter is intimately related to the former and is known to be particularly 
important in the threshold region. In the bW^bW~ production, there are non-resonant 
diagrams where bW^ and bW~ are not produced from the decay of t and t. Often the non- 
resonant diagrams are omitted in the studies of ti productions, since they are suppressed in 
the events where both t and i are nearly on-shell. In the threshold region, however, either 



of t and t tends to be off-shell due to restricted phase-space and the binding effects ||2^ , and 



the non-resonant diagrams can give non-negligible contributions compared to the resonant 



{ti) diagrams [31|. Since these contributions interfere with each other, all the diagrams have 
to be taken into account at the amplitude level. Moreover, because of the requirement by 
unitarity, we also have to ensure a consistent treatment of the finite decay width of the 
top quark in our framework. We will discuss these points within our framework, which 
includes the bound-state effects as well as the non-resonant diagrams, in connection with 
a known problem regarding a gauge cancellation. 

In order to compute numerically various kinematical distributions at hadron colliders, 
we develop a MC event-generator, which is adapted to the MadEvent [^, ^] environment ^^ . 
We include the bound-state effects in the hard-scattering part of the LO event-generator, on 
the basis of our theoretical framework. The ISR and/or final-state radiation (FSR), which 
are of importance at hadron colliders, are incorporated via the parton-shower approach. 
Since the parton shower does not alter the normalizations of the cross sections at the 



partonic level, we will complement the overall normalizations, known up to NLO |2C, 21 1, 
by multiplying the cross sections in the individual channels with the so-called "i^-factors." 
We note, however, that our aim here is to construct a generator valid only up to LO 
with respect to the differential distributions, in this first attempt to include the bound- 
state effects. Compare with the existing NLO event-generators, such as MCSNLO [^, ^ 
and POWHEG [^], which realize a consistent treatment of perturbative corrections for any 
process and any phase-space point. 

Using the generated events, we study the bound-state effects on the top quark differ- 
ential distributions at the LHC, focusing on the events in a relatively low m^j- region. A ti 
pair gains a binding energy due to exchange of Coulomb gluons between them. This effect 
tends either of t and t to be off-shell below the threshold, and the effect remains even a few 
tens GeV above the threshold, due to the large width of the top quark. We will quantify 
this picture through detailed examinations of the top quark differential distributions. 

The paper is organized as follows. In Sec. 2, we give a theoretical framework for 
computing the amplitudes for top-quark pair-production at hadron colliders, incorporating 
the bound-state effects (Sec. 2.1), the finite width effects (Sec. 2.2), and the ISR effects 
and X-factors (Sec. 2.3). In Sec. 3, we present numerical studies for various kinematical 



^The Fortran code for the event generator including the bound-state corrections is available at 
http : //madgraph . kek . jp/'yokoya/TopBS/. 
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Figure 1: Feynman diagrams for gg -^ tt and qq — >■ tt at the tree- level. 

distributions in tt production, using the MC simulation which implements the ingredients 
explained in the previous section. In Sec. 4, we summarize our results. To avoid complexity 
in the main body of the paper, several detailed discussions are presented in the Appendices. 
In App. A, we identify the ii Green function in a Feynman amplitude. In App. B, we derive 
the off-shell suppression factor. In App. C, the color decomposition of the amplitude is 
explained. In App. D, we examine the leptonic decays of W^s from top quarks with and 
without spin correlations. 

2. Inclusion of Bound-state Effects 

In this section we present a theoretical investigation of how to include the ti bound-state 
effects in the matrix elements for gg — t- bW^bW" and qq — t- bW'^bW" . In particular we 
include the effects such that the amplitude is correct in the leading-order approximation 
both in the tt threshold region and in the high energy region. Inclusion of several different 
effects is explained in steps: In Sec. 2.1 we explain how to incorporate the bound-state 
effects; in Sec. 2.2 important higher-order effects of the large top-quark decay-width are 
incorporated; in these subsections, we consider only the partonic 5-matrix elements. In 
Sec. 2.3 we incorporate the ISR effects and the /f-factors in the corresponding partonic 
differential cross sections. 

For later convenience, we divide each amplitude into two parts, the tt (double-resonant) 
part and the non-resonant part, as 

M^^\l -^ bW+bW) = M^^{I -^ tt^ bW+bW-)+M^^){I -^ bW+bW-) , (2.1) 

where I = gg or qq represents initial-state partons, and c represents the color (c = 1 and 
8 for the singlet and octet, respectively) of /, or equivalently, of bb in the final-state. The 
first term on the right-hand side represents the sum of the diagrams which contain both 
t and i as an intermediate state. This part of the amplitude consists of I — )• tt processes 
followed by subsequent decays of t and t. The second term represents the sum of the rest of 
the diagrams, which consists of single(-top)-resonant diagrams and non-resonant diagrams. 
Fig. H] shows the tree-level Feynman diagrams for the processes gg — t- tt and qq — )• tt. Some 
examples of the tree-level diagrams included in each part for I = gg are shown in Fig. |2|. 
In general each part is gauge-dependent. In this paper, we work in Feynman gauge for 
SU{3)c and in unitary gauge for the broken electroweak symmetry. 

In computing the tree-level Feynman diagrams which contain the top-quark propaga- 
tors, we include the (on-shell) top-quark decay- width F^ in the propagator denominator 
as 

■SFiPt) = -^ ^—. — w- ■ (2-2) 

pf -mf + initTt 






Figure 2: Typical Feynman diagrams for the double-resonant (left), single-resonant (middle) and 
non-resonant (right) contributions in the gg — ^ bW^bW~ process. They belong to -Mj^ , -Mn/ and 
A^nr , respectively. The term "resonant" is used to refer to the t or t quark propagator (shown with 
double line) which can become close to on-shell. The blob in the left diagram represents the first 
three diagrams in Fig. |l| 

2.1 LO cross section valid from threshold to high energies 

In this subsection we include the bound-state effects in the tt amplitude -A4Jf • We consider 
the narrow-width limit of top-quark width in this subsection. Namely, we take into account 
only the leading contributions as Tt — )• 0. Important subleading effects by the finite top- 
quark width will be investigated separately in the next subsection. 

We start by reviewing the conventional method for including the bound-state effects 
in the matrix element (or in the fully differential cross section) of e"*"e~ — t- tt — t- bW^bW~ 
close to the threshold of tt pair productions. At the leading-order, this is achieved by 
multiplying the tree-level amplitude corresponding to the diagrams e~^e~ — )• tt — t- bW~^bW~ 

lit 



by an enhancement factor as |25, 



tree 



M{e+e- ^tt^ bW^bW") = M{e+e- ^ tt ^ bW+bW 

"" Go{E + iTt,p) 
Here, the non-relativistic Green functions are defined by^ 

iE + ^Tt)-i-^ + V^^i^{r) 



(2.3) 



mt 



G^''\E + iTt,f^ = 5^{f^, (2.4) 

G^''\E + iTt,p)= I d^re-'^-'^&''\E + iTt,r). (2.5) 

E = \fs—2mi is the cm. energy measured from the threshold; rrij is the pole mass of the top 
quark; r denotes the relative coordinate of t and t, while p denotes the three-momentum of 
t (or minus the three- momentum of t), both defined in the cm. frame; VQ^j-)(r) is the QCD 
potential between the tt pair in the color-singlet (c = 1) or color-octet (c = 8) channel. In 
e+e~ collisions, tt pairs are produced in the color-singlet channel, hence c = 1 in Eq. (p.3|). 
The free non-relativistic Green function Go{E + iTt,p) is obtained from G^'^'{E + iTt,'p) by 



•''There are two tree-level diagrams for e+e~ -^ ti ^ bW^bW~ with 7 and Z boson intermediate states. 
M{e^e~ -^ ti ^>- bW^bW~)trcc denotes the sum of them. 

^In this study, G^'^\E 4- iTt,p) is computed numerically by solving the Schrodinger equation in coor- 
dinate space and taking Fourier transform E3]. Alternatively, one may solve the Schrodinger equation in 
momentum space directly |26[. 



ric) 



setting Vqqy){^) — ^ 0. Formally the above Green function can be expressed as 



G(^)(s + iri,p) 



p 



E + iTt-pymt-V^^^{r) 



(2.6) 



using an operator notation in quantum mechanics. By definition, the above Green function 
contains only the 5"- wave contributions. 

There are two methods to compute the total cross section for e~^e~ — )• tt — )• bW~^bW~ 
incorporating the bound-state effects. One way is to integrate the absolute square of the 
matrix element given in Eq. (2.3) over the phase-space of the final bW~^bW~ state. The 
other method is to use the optical theorem (unitarity relation) and take the imaginary 
part of a current-current correlator (two-point function). At the leading-order, the latter 
method leads to the formula l27^, PStl: 



(Ttot(e+e- ^tt^ bW^bW 



o-tot(e^e 



Im 



^'■jtree 

&^\E + iTt,r = 0) 



Im 



Go{E,r = 6) 



(2.7) 



where the Green function in coordinate space is defined in Eq. (p.4|). o"tot(e'*'e~ — t- tt)tree 
denotes the Born cross-section for the production of on-shell top quarks. Note that in the 
denominator we set Fj to zero in Gq, whereas F^ is retained in the denominator in Eq. (p.3|). 
The different treatment of F^ is because in Eq. (p.3p we use the tree-level amplitude with 



unstable top quarks (in the intermediate state), whereas in Eq. (2^) we use the tree-level 
cross section of the on-shell top quarks (in the final state). 

Both formulas ( ^.3| ) and (2^) incorporate all the leading-order corrections ~ (as/f^)"' 
in the threshold region E <^ rrit. On the other hand, the formulas are not valid at higher 
cm. energies E ^ mt, since relativistic corrections ~ /?", which grow with energy, are 
neglected in these formulas. 

Now we turn to the partonic cross sections for the top-quark productions in hadron 
collisions, / — )• tt — )• bW~^bW~ with I = gg and qq. Unlike e^e~ collisions, the collision 
energy of the initial state cannot be fixed. Thus, we need to consider both threshold and 
high-energy regions. The formulas which we propose, valid in both regions within the 
leading-order approximation, can be summarized as follows: 



(1) The tt amplitude for / 

M^Sil^ti 

with 



tt 



bW'^bW is given by 



bW+bW- 



M^H 



E' = E + 



E' 
Arrii 



tt 



bW+bW-)tree 

G'^^XE' + iVt^p) 
Go{E' + iTt,p) ' 



(2. 



(2.9) 



ic)i 



tt 



The Feynman diagrams which contribute to the tree-level amplitude M,t {I 
hW^hW~)i^ee are those shown in Fig. |^, after attaching the decay vertices t — )• bW'^ 
and i — >• hW~ to each diagram. There are both color-singlet (c = 1) and color-octet 
(c = 8) channels in the case I = gg, while there is only the color-octet channel in the 



case / 



Here, E is defined from the tt invariant-mass tti^j as E 



771,4 



The only essential difference from the corresponding formula for the e~^e 
Eq. (P^), is the use of the modified energy E' [Eq. (|2.9D ] instead of E. 



„„ 2777t. 

collision, 



(2) We may compute the it invariant-mass distribution by integrating the absolute square 
of the above amplitude \M]^ {I — )• tt — )• bW'^bW~)\'^ over the bW~^bW~ phase-space 
for each fixed mff. Alternatively we may obtain the formula for the ti invariant- 
mass distribution using the optical theorem, similarly to the e"'"e~ — )■ ti case. In the 
leading-order approximation, this reads 



<tS(/ ^tt^ bW^bW-) = a. 



Im 


'&^ 


^E' + iTt,f = 


0) 


totl-' ' "J tree '^ 


Im 


Go{E',r = 0) 





As in the e~^e 



(2.10) 

ti case, d'l^l{I — )• tt)ti.ee is the Born cross-section for the on-shell 
top quarks; accordingly Tt is set to zero in GQ{E',f = O). We note that the ti 
invariant-mass distribution obtained from Eq. ( ^.10[) does not exactly coincide with 
that obtained by integrating \Mlf{I -^ ti ^ bW^bW~)\'^ over the bW~^bW~ phase- 
space; the difference is 0{rt/mt) and will be discussed in the next subsection. 

In the following we sketch our theoretical consideration which led to the above formulas 
( |2.8| ) and ( p. 10 ), where some details are relegated to App. ^. As explained in that appendix, 
part of the Feynman amplitude for / — t- tt — )■ bW^bW^ can be identified with a Green 
function that dictates the time evolution of the tt system. (Such an identification is possible 
in all kinematical regions.) In the cm. frame of ti, for the initial-state | i ) and final-state 
I / ) of the ti system, this Green function can be written formally as " 

^f\ ^TX^I*)' (2.11) 

where rriit is the cm. energy of tt (tt invariant-mass). The full QCD Hamiltonian is denoted 
by H. Because of this property, the amplitude for I ^ ti ^ bW^bW~, which incorporates 
the ti bound-state effects, can be obtained from the tree-level amplitude by multiplying an 
enhancement factor: 



Mlyil ^tt^bW^bW- 



M^fil^ti 



bW+bW-)tree 

1 



mti - H + ir 



P 



1 



ruu - Hq + iVt 



(2.12) 



^As already stated, treatment of the top-quark width is correct only in the leading order. 



Hq denotes the Hamiltonian H after setting as — )■ 0, i.e. the free Hamiltonian. As before, 
r denotes the relative coordinate of t and t, while p denotes the three-momentum of t, 
both defined in the cm. frame of the ti system. Corrections to Eq. (|2.12| ), which come 
from the non-resonant part of the amplitude (i.e. that vanish as Tt — )■ 0), are neglected; 
see Appendix ^. 

In the above equation, we have taken advantage of the fact that we work in the leading- 
order approximation and set the initial-state of the Green functions as |r = 0). This follows 
from the following consideration. Naively, the ti system cannot be regarded as being created 
by contact interaction (i.e. at the same point r = 0) in the t- and it-channel diagrams of 
gg —7- tt, in which the top-quark is exchanged between the initial gg (see Fig. |l|) . In the non- 
relativistic region, however, we may set p — t- in the t- or u-channel top-quark propagator 
in the leading-order approximation. The denominator of the propagator effectively reduces 
to a constant close to the threshold, so that tt can be regarded as being created at the same 
point. On the other hand, in the relativistic region the enhancement factor in Eq. ( 2.12| ) 



reduces to 1 + 0{as)- Hence, it is justified to evaluate the Green functions with the initial- 
state \f = 0) within our present approximation. 

Since our aim is to include the leading-order contributions both in the relativistic 
and non-relativistic regions, the full form of the Hamiltonian is not necessary. In the 
region where t and t are relativistic, rriu; — 2mt ^ rrit, the leading-order contribution in the 
Hamiltonian reads 



H = 2^p2 + ^2 ^ 0(^^) ^ (2.13) 

as shown in Appendix ^. It is nothing but the sum of the energies of free on-shell t 
and t. The above equation also indicates how the next-to-leading order effects enter the 
Hamiltonian. On the other hand, in the non-relativistic region, E = mti — 2mt <^ mt, the 
leading-order contributions in the Hamiltonian can be written explicitly as 



H = 2mt + 



I + ^cbM 



l + 0(a„/3) . (2.14) 



It is indicated that the next-to- leading order corrections enter as 0{I3) relativistic correc- 
tions or 0{as) corrections. 

A natural choice of the Hamiltonian, which incorporates the leading-order contribu- 
tions in both regions and smoothly interpolates these regions, is given by 



i^LO = 2Jp^ + ml + Vi'l^ir). (2.15) 



In fact, it is well known that, when computing higher-order corrections in Coulombic 
bound-state problems, part of them (relativistic corrections) contribute exactly in the above 
form [37]. Thus, in principle, one may determine the enhancement factor in Eq. ( p. 12 ) using 



the above Hamiltonian. 

Due to technical reasons, however, we use an alternative form of the enhancement 
factor, which is equivalent within the present approximation. By substituting E = rriii — 



2mt to the on-shell relation 



nitt = '2\/p + m 



"i' 



one finds that 
Therefore, if we define 



E + 



E^_p^ 
Amt rrit 



E' + iVt-{p-^/mt + vi^l^{r)) 



f =0 



(2.16) 



(2.17) 



(2.18) 



with E' defined by Eq. (2.9), the position of the pole oiG^'^'{E' + 1^1,"^) is the same as that 
of {p\{m^i — Hj^o + iVi)~'^\f = 0) in the limit ag — >• 0. Furthermore, in the non-relativistic 
region, evidently G^'^\E',p) is the same as (p|(wijf — -^lo + ^r*)"^!^ = 0) in the leading- 
order approximation. Hence, we may compute the matrix element by the formula Eq. (p.^). 
One may be worried that, although the replacement E ^ E' correctly accounts for the 
pole position, the residue of the pole may be altered significantly from that of Eq. ( 2.12| ) in 
the relativistic region. This is not the case, since the change of the residue will be canceled 
in the ratios of the Green functions in Eqs. i \i.8\j and ( 2.10| ). 

The advantage of using G^^'{E' A-iTt^p) is that one can obtain it from the conventional 
non-relativistic Green function with a minimal modification E ^ E' = E + E"^ /{Anit). In 
particular, properties of G^'^' are fairly well known. 

Let us comment on the dependence of the Green function on the top-quark width Tt in 
Eq. ( p. 18 ). In Eq. (|2.2|) the shift of the pole position of the top-quark propagator due to the 



finite top-quark width can be incorporated simply by a replacement mj — )• mi — irritTt. If 
we apply it to Eq. ( p.l6| ) , one finds that iVt will be added to the left-hand side of Eq. ( 2.17| ) . 
Hence, inclusion of Tt as in Eq. ( |2.18| ) is correct in the leading-order approximation. 

Throughout our analysis, we include an important subleading correction to the bound- 
state effects, in order to make our analysis more realistic. This is the NLO (1-loop) cor- 
rection to the static potentials ^^^^(r) between the tt pair. The NLO potential reads |^8| 



K 



(c) 



QCD 



(r;^B) = C(^) 



as{^^l 



1 + 



«s(^E 
A-n 



{2/3o[ln(^Br)+7£] + 4'^} 



(2.19) 



with 



C(i) = -Cf, C(8) 


Ca ^ 
- 2 Cp, 






(2.20) 


/3o = ^Ca - -n„ 


af = af - 


3V 


10 


(2.21) 



for the MS coupling. Here, ^e = 0.5772... denotes the Euler constant; Cp = 4/3 and 
Ca = 3 are color factors. The QCD potential is renormalization-group invariant, and 
we evaluate the above expression at the Bohr scale of /iB = 20 GeV and with n„ = 5 
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Figure 3: Partonic tt invariant- mass distributions for gg 
channel. The green solid line is calculated with E' — E 
dot-dashed line is calculated with E instead of E' . The black dashed line represents the Born cross- 
section. The dotted lines are those for the gg — > tt process (production of on-shell top quarks). 



(Q,(/is) = 0.153). 



Now we perform a few tests of our formulas, Eqs. ( |2.8D and ( 2.10| ). First we examine the 
impact of the replacement E ^ E' . The tt invariant-mass distributions da/drritt computed 
with these formulas are compared with da/dm^f computed by the formulas valid only in 
the threshold region, namely Eqs. ( |2.8| ) and ( p. 101 ) after we replace E' by E. In Fig. ^ 
we plot a for gg — t- ti in the color-singlet channel." The green solid and dotted lines 



are those computed using Eqs. (2.8) and (2.1C), respectively, while the red dot-dashed 
and dotted lines are those computed with the same formulas but after the replacement 
E' — )• E. For comparison, the Born cross-sections [using Eqs. (2^) and ( |2.10| ) but without 
the enhancement factors] are also plotted with the dashed and dotted black lines. All 
the cross sections are computed with rrit = 173 GeV and Tt = 1.49 GeV (the tree-level 
top-quark decay-width) . 

The difference between the ti invariant-mass distributions using Eqs. ( p.8[ ) and ( 2.1C| ) 
(solid and dotted green lines) is due to 0{Tt/mt) corrections. The replacement E' ^ E 
in Eq. 



.101) changes the tt invariant-mass distribution slightly above the tt threshold; 
compare the green dotted and red dotted lines. The difference between the two cross 
sections is about 2.5% in the large m^j region. 

On the other hand, the effect of the replacement E' ^ E in Eq. ( p. 8] ) is much more 
pronounced above the ti threshold. There exist a large enhancement which amounts to 



"in Sees. 2.1 and 2.2, the partonic ti invariant-mass distribution (before including the effects of ISR 
and parton distribution function) is proportional to the partonic total cross-section and delta function, 
da/dnitt oc a5{s — m.%)\ c.f. eq. (17) of |20|. Hence, we plot (t(s = m^) instead of da/drriti in these 
subsections. 
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Figure 4: Ratios of the partonic cross sections for gg — > hW^hW^ and gg -^ tt, in the color- 
singlet channel at rTijj — 360 GeV. These are plotted as functions of Ft- Green solid, red dashed 
and black dotted lines are those with the modified energy E', with E instead of E' , and for the 
Born cross-sections, respectively. 



conditions satisfied by the pole positions of the t and t propagators contained in A^Jr! 
by the pole position contained in Go {E + iVt , p) ■ [Note that al^^. ^.^^^ 



nearly a factor of two around yn^f = 400 GeV; compare the green solid and red dot-dashed 
lines. The origin of this large enhancement can be identified with a mismatch of the on-shell 

(^tot.tree i^ Eq. ( ^ ) docs not 
contain the t or i propagator, so that this mismatch problem does not occur when we replace 
E' by E in Eq. ( 2.1C| ); compare the green and red dotted lines.] In fact, the mechanism of 
this abnormally large deviation is closely tied to a characteristic bound-state effect on the 
invariant-mass distributions of the bW^ and bW~ systems. We will investigate this issue 
in detail in Sec. 3, in which we examine closely the differential distributions. Nevertheless, 
even without going into these details, the present comparison clearly shows the necessity of 
a proper treatment of the relativistic kinematics, when we include the bound-state effects 
to the fully differential cross section of the process / — )• tt — )• bW^bW~. 

Since our formulas are correct in the narrow width limit Tt — )• 0, the unitarity relation 
should be restored in this limit. In order to check this, in Fig. ^ we plot the ratios of 
da /drriii computed using Eqs. (^) and ( p.lOj ) as we vary the value of Tt at a fixed nitt of 
360 GeV.** We confirm that the ratio approaches to unity as Tt is reduced, in the case that 
we use our relativistic formulas (green solid line) or in the case that we use the tree-level 
cross sections (black dotted line). In sharp contrast, the ratio does not approach to unity 
in the case that we replace E' by E (red dashed line) due to the mismatch problem. It 
shows invalidity of the non-relativistic approximation far above the threshold, especially 
for the fully differential cross section. 



**As we vary the value of Tt in the t propagator and the Green functions, the value of the weak gauge 
coupling constant gw in the thW vertex is varied consistently, such that the tree-level top-quark width 
takes the correct value. 
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Figure 5: Partonic if invariant- mass distributions for gg — > tt in the color-singlet (red), color-octet 
(blue) channels, and the sum of them (black). The solid lines are calculated by multiplying the 
bound-state enhancement factor to the whole ti amplitude. The dashed lines are calculated as the 
sum of contributions from all J's, where only the J = amplitude is multiplied by the enhancement 
factor. The dotted lines represent the Born cross-sections corresponding to the above lines. 



As is well known, the leading-order bound-state effects in the ti threshold region are 
contained in the S-wave part of the amplitude. In the case of gg — )• tt, the 5-wave contri- 
butions reside in the J = amplitude both for the color-singlet and color-octet channels. 
Hence, it may be more appropriate to include the bound-state effects only in the J = 
amplitude, rather than multiplying the whole ti amplitude by the enhancement factor as 
in Eq. (2.8). Theoretically, the difference between the two prescriptions is subleading. We 
examine this feature by comparing the ti invariant-mass distributions computed in both 
ways. In Fig. |5|, we plot the ti invariant-mass distributions for gg — t- ti process. Each solid 
line represents the cross section computed using Eq. ( p.lOD , namely, the whole Born cross- 
section (the sum of the Born cross-sections for all J's) is multiplied by the enhancement 
factor. Each dashed line represents the sum of the cross sections for all J's, where only the 
J = cross section is multiplied by the enhancement factor. The red (solid and dashed) 
lines represent the cross sections for gg — )■ ti in the color-singlet channel, while the blue 
(solid and dashed) line represent those in the color-octet channel. 

The cross section in the singlet channel is more enhanced if we use the overall prescrip- 
tion, Eq. ( p. 10 ), since the force between the color-singlet ti pair is attractive and hence 
the enhancement factor is larger than one. The difference of the two prescriptions is siz- 
able only above the ti threshold and becomes maximal around ttt-^j ~ 400 GeV, where the 
difference is about 7%. On the other hand, the cross section in the octet channel is more 
reduced if we use the overall prescription, since the force between the color-octet ti pair 
is repulsive and the enhancement factor is (slightly) less than one. The difference of the 
two prescriptions is at most 2%. The black lines (solid and dashed) represent the sum 
of the cross sections for gg — t- ti in the above two channels. The difference of the two 
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prescriptions in this case is at most about 1%, since the differences have opposite signs in 
the two channels and are largely canceled. Thus, the difference of the two prescriptions 
is rather small and much smaller than other subleading corrections which we neglect in 
our analysis. Furthermore, we have checked that the above tendencies are not changed 
significantly by the ISR effects. Therefore, for simplicity of our analysis, we will adopt the 
overall prescription in the following analysis, namely, we will not decompose the amplitude 
into different J's. 

In the case of qq — t- it, there is only the J = 1 color-octet channel at tree level. Hence, 
the enhancement factor multiplies the whole amplitude also in this case. 

2.2 Effects of large Vt 

In this subsection we describe how we incorporate part of the subleading corrections that 
are induced by the large top-quark width. As an inevitable consequence of numerically 
integrating the fully differential cross sections for / — t- it — >• bW^bW~ , there is a significant 
phase-space-suppression effect. We partly compensate this effect, which is related to a 
gauge cancellation inherent in the inclusive cross section. 

First, we briefly review existing theoretical studies on the treatment of the top-quark 
width, in the cases with and without bound-state effects. In the latter case, many schemes 
have been proposed for incorporating the top-quark width. The use of the top-quark prop- 
agator in Eq. ( |2.2| ) is called the fixed-width scheme (FWS). It is widely used in simple 
analysis of the cross sections which include the top quark as an unstable intermediate par- 
ticle. It is known, however, that subleading electroweak effects are not properly treated in 
this scheme. At present, the complex-mass scheme (CMS) [^9[ seems to be most advanced 
from a practical point of view, due to the simplicity of its implementation. In fact, for 
the process e~^e~ — )• W^W~ , which is kinematically similar to ti productions, the fully 
differential cross section has been computed incorporating the effects of M/^-boson width 
with NLO accuracy in this scheme, basically in all kinematical regions ^^. For tt pro- 
ductions in hadron collisions, the fully differential cross sections is computed incorporating 
top-quark width with LO accuracy in CMS, and various differential cross sections in differ- 
ent schemes were compared |30|. In particular, the study has shown an agreement within 



errors between all the calculated cross sections in CMS and in FWS. (FWS is simpler but 
less sophisticated than CMS.) 

Regarding ti productions in the threshold region including the bound-state effects, 
studies on the finite- width effects are most advanced in the total cross section for e^e~ — )• 
ti. The finite- width effects have been incorporated with NNLO accuracy* [O, H2l, using 



the velocity-Non-Relativistic QCD (vNRQCD) effective field theory framework |43, 44|. 



Recently an NLO correction to the total cross section arising from the single-top resonance 



region has been pointed out and computed in |31|, using unstable-particle effective field 



theory |45, 46 1. On the other hand, in the corresponding fully differential cross section the 



width effects are incorporated only up to NLO accuracy [47| (apart from the contributions 



from the single-top resonance region). MC generators, developed specifically for simulation 



"In the tt threshold region, it is customary to count Tt/nit ~ 0{aw) ~ 0{a^). 
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studies in the threshold region of the e^e — )• tt process, have incorporated both bound- 



state effects and finite-width effects in the LO approximation [ 22 , ^] . For ti productions in 
hadron cohisions, only the ti invariant-mass distributions have been computed with NLO 
accuracy, incorporating the bound-state effects, finite-width effects and ISR effects in the 
threshold region [20, p| ]. 



One effect is known to be particularly important in computing the fully differential 
cross sections in the threshold region of ti productions. It is the phase-space-suppression 



effect 1 25, 26, H3|, which is formally an NNLO effect of the top-quark width, but it seriously 



modifies the shape of the sharply rising S'-wave cross section as a function of mtj, after 



integrating the differential cross section over the final-state phase-space |25]. Let us briefiy 
explain this effect. The ti cross section starts to rise below the ti threshold rriii = 2mt 
as a result of formation of ti resonances. This means that the dominant kinematical con- 
figuration is such that one of t and i is on-shell and the other is off-shell. Therefore, the 
phase-space of bW which decayed from the off-shell t or t is reduced as compared to the 
on-shell case. This suppresses the production cross section, and this effect is automatically 
incorporated if we integrate the LO differential cross section numerically over the phase- 
space of the final bW~^bW~ . A remarkable feature is that there is another effect at NNLO 
which exactly cancels the phase-space-suppression effect, for the integrated cross section 
at each ti invariant-mass [^, Q . This is the Coulomb-enhancement effect due to gluon 
exchanges between t and b (decayed from i) and between i and b. The cancellation is guar- 
anteed by gauge invariance (Ward identity) . It protects the ti resonance widths from being 
determined by gauge-dependent off-shell width of the top quark. Consequently the only 
surviving NNLO effect to the ti resonance widths turns out to be the time-dilatation effect 
due to the relative motion of t and i inside the resonances, which is gauge independent. t 

Thus, we face a problem when we compute differential cross sections in the tt threshold 
region by a MC generator: The phase-space-suppression effect is automatically incorpo- 
rated, while the Coulomb-enhancement effect due to gluon exchanges between t and b (or t 
and b) is difficult to incorporate in a MC generator. (This is not yet achieved even in theo- 
retical computations of the e^e~ — )• ti differential cross sections.) Our prescription in this 
study is only effective. Since we know that the phase-space-suppression effect is canceled 

— — (c) — J_~ 

in the inclusive tt cross section, we multiply the tt amplitude ^A^^ {I -^ tt ^ bW^bW~) 
by an enhancement factor such that the phase-space-suppression factor is canceled. In 
addition, we include the non-resonant diagrams, which are formally 0{Tt/mt) compared 
to the LO amplitude. 

Hence, our full amplitude at the parton level is given by 



-^it+SH/- (^ ^ bW+bW) = m'^ + M^J (2.22) 



^In principle, the momenta of t and t can be determined from the final state, in the LO approximation. 
Hence, the relative motion is a gauge-independent quantity. 
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(2.23) 
(2.24) 



(c) 

where A^)r is defined by Eqs. 



Ac) 



and (12.91), and TW^r denotes the sum of the tree- level 



non-resonant diagrams. The factor in the square bracket represents the inverse of the 
phase-space-suppression factor (see App. B for the derivation), where Tt{st) denotes the 
running top-quark width in unitary gauge evaluated at the top-quark invariant-mass Sj; its 



Ac) 



Ac) 



explicit form is given in App. B, Eq. (B^). At large tt invariant-masses, the bound-state 



effects diminish, namely M..T — )• A^lf i > hence the above bW~^bW amplitude is defined 



Ac) 






such that it reduces to the tree-level amplitude, ^A, _f_T 

At the differential level, the above treatment of the cancellation of the phase-space- 
suppression effect is only effective, since the Coulomb-enhancement effect does not cancel 
the phase-space-suppression effect at each kinematical point. Nevertheless, we consider 
that a higher priority should be given to the gauge cancellation mechanism that is inherent 
in the inclusive cross section. We also note that the replacement E —^ E' in the Green 
function, Eqs. ( |2.8D and (|2.9|), automatically incorporates the time-dilatation effects to the 
resonance widths. "^ 

In Fig. P, we compare the tt invariant-mass distributions for gg — t- tt and qq — t- tt in 
all the channels, which are computed by integrating the following four cross sections over 
the bW^bW~ phase-space, where the tt invariant-mass rritt is defined as the invariant-mass 
of the final bWHW system: (i) |Xjf P (black dotted), (fi) |A^Jf P (blue dashed), (iii) 
|A^jj + -A4nr P (^6*^ dot-dashed), and (iv) the absolute square of our formula Eq. ( |2.22| ) 
(green solid). Comparing the distributions for (i) and (ii), or, (iii) and (iv), we see that the 
phase-space-suppression effects are sizable especially close to the threshold of gg — )• tt in the 
color-singlet channel. This is consistent with the explanation given above. In particular, in 
each figure the difference between (iii) and (iv) is hardly visible at large m^f. We may also 
compare the distributions for (i) and (iii), or, (ii) and (iv), to see the contributions of the 
non-resonant amplitude. The contributions become comparatively larger at high energies 

(c) 

for gluon-fusion channels, since the contribution of the s-channel diagram in M)t decreases, 

(c) 

while contributions of the single-resonant diagrams in A^V increase. On the other hand, 
for qq channel, contributions from the non-resonant diagrams are small everywhere. 

In Fig. I^a) we plot the tt invariant-mass distributions for gg — )• bW^bW~ (sum of 
the color-singlet and octet channels), using Eq. ( p. 22 ) after integrating over the bW~^bW~ 
phase-space (red solid line) . The tt invariant-mass distribution at the Born level for gg — )• 



''The relation E' + iVt = p^/mt (corresponding to 



m. 



m. 



imtFt) is relativistically correct, so that 



^2 



the time-dilatation effect enters the lifetime of the tt system. It can also be seen by the fact E' + iTt — -2— ~ 
E + iFt-^ + ^i^- iVt] ~ E + irt(l- ^) - ^ + -^, where (l - ^) represents the 

-i 

time-dilatation effect and -r- 

4m: 



iFtj ~ E + irt(^i 

represents a relativistic correction 
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Figure 6: Partonic tt invariant-mass distributions for gg — )■ tt and qq — > tt in the individual 



channels. Four lines in each figure are computed from (i) jA^^^ | (black dotted), (ii) [A^^^ | (blue 
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iic) 



<(c)|2 



(c) 



<■(==) 1 2 



dashed), (iii) \Ml^ + A^nr | (red dot-dashed), and (iv) {Mj.^ -I- A^nr | (green solid). The figures 
in the right show the magnification of the threshold region. 
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Figure 7: tt invariant-mass distributions for the processes gg -^ hW^bW^ (color-summed) and 
qq — > bW^bW~. The red solid line is plotted using Eq. ( ^.22 ) for the partonic amplitude, while 
the black dashed line represents the Born-level cross section. 



bW^bW^ in FWS is also plotted (black dashed line). The difference of the two lines signifies 
the bound-state effects, after including the non-resonant diagrams and compensating the 
phase-space-suppression effect. The enhancement of the distribution by the bound-state 
effect is visible not only in the threshold region but also up to about ?7ijf = 500 GeV. The 
enhancement factors are about 1.05 and 1.02 at rrift = 400 GeV and 500 GeV, respectively. 
The tt invariant-mass distributions for qq — t- bW^bW~ are also plotted in Fig. 0(b). The 
enhancement factor is smaller than unity because of the repulsive force, whose values are 
about 0.96 and 0.98 at m^f = 400 GeV and 500 GeV, respectively. 
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Figure 8: rrift dependence of the iiT- factors for tt production in gg color-singlet (top), color-octet 
(middle) and qq (bottom) channels at the LHC y/s — 14 TeV. The dashed, solid and dotted lines 
are obtained with /i/; = fip = unit with k = 0.5, 1 and 2, respectively. 

2.3 Inclusion of ISR effects and ET-factors 

In this subsection we explain how we incorporate the ISR effects in the cross sections in 
our framework. In addition, we determine the i('-factors to match our predictions for the 
tt invariant-mass distributions to the available NLO predictions. 

In hadron collisions, it is important to include the ISR effects. In our framework, 
they are incorporated by connecting the differential cross sections computed from the 



matrix elements Eq. (2^) to a parton-shower simulator such as PYTHIA |51] or HERWIG [p2| . 
In addition, we include "X-factors" as the normalization constants of the cross sections 
for / — )• bW^bW~ in the individual channels.^ The i^-factors are determined such that 
the ti invariant-mass distribution for each channel in the threshold region matches the 
corresponding NLO prediction in the threshold region. We also extrapolate these K- 
factors to the large m^f region. The main reason to do so is a lack of the NLO predictions 
in the large mfj- region for the individual channels: The present theoretical prediction for 
the NLO ti invariant-mass distribution in the large m^j region is provided numerically only 
for the color-summed cross section for on-shell top-quark productions |l^. Theoretically, 
by naive extrapolation of the -ftT-factors, we reproduce the double-logarithmic terms of 
the cross sections correctly in the large m^f region, due to the universal structure of soft- 
gluon emissions; on the other hand, we do not reproduce the single-logarithmic and non- 
logarithmic terms. 

The NLO corrections to the ti invariant-mass distributions in the threshold region 
are known for the individual channels [^, ^]. The corrections are given in terms of the 
hard-correction factors and the gluon radiation functions. The major difference of the 
predictions of [20| and |21| is that in the latter predictions contributions from high yl 

^Note that the parton shower simulators incorporate ISR effects by way of stochastic processes, and that 
the ti invariant-mass distributions handed to the simulators at the parton level are not affected by the ISR 
effects. 
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(the cm. collision energy of the initial partons) are included more accurately. " Hence, we 
use the latter predictions to compute the iT- factors". We can determine the i^- factors by 
taking the ratios of these NLO partonic cross sections and our (LO) partonic cross sections 
given in Sec. 2.2. Since the NLO cross sections |2C, 21 1 do not include contributions from 



(c) 

non-resonant diagrams A^nr , accordingly we incorporate only the contributions from the 
resonant diagrams A^)f in the LO cross sections when we compute the i^-factors. (In most 

— (c) 

of the threshold region, the effect of A^V is irrelevant in any case, since resonant diagrams 



dominate.) Furthermore, in calculating the i^-factors, we use the CTEQ6M PDFs [56] and 
the 2-loop running of the strong coupling constant a^ for the NLO m^f distribution, while 
the CTEQ6L1 PDFs and the 1-loop running of q^ are used for the LO distribution. (We 
find that the i^-factors obtained by using the MSTW2008 PDFs [ p7| ] are quite similar.) 

In general, the i^-factors depend on m-jf. We first examine ?TT,^j-dependences of the 
ii'- factors as we choose different renormalization and factorization scales, fiji and ^ip- The 
renormalization scale /iR enters the NLO formula as the scale of the strong coupling con- 
stant and also through the logarithmic term in the hard-vertex function; see Eq. (3.2) of 



1 21]. On the other hand, the factorization scale ^p enters the NLO formula as the scale of 
the parton distribution functions (PDFs) and through the terms with In ijnl^/ ^i^p) in the 
gluon radiation functions; see Eqs. (3.4-3.7) of [^]. We find that, the TTi^j-dependences of 
the i^-factors can be relatively flat in the threshold region, with appropriate choices of jjlr 
and ^,F- In this case, extrapolation of the iC- factors from the threshold region to the high 
mil region can be performed trivially. Indeed, for simplicity of our analysis, we take the 
-RT- factors to be independent of m^f . In Fig. |^, we plot the i^- factors of the ti invariant-mass 
distributions in the individual channels at the LHC with ^/s = 14 TeV, in the cases that we 
choose the scales as hr = fip = nrnt with k = 0.5, 1, 2. As can be seen, the m^j-dependence 
of the iiT- factors are mild. We have also examined the i^-factors corresponding to the LHC 
with -y/i = 7 TeV and Tevatron with ^/s = 1.96 TeV; we find that the X-factors are only 
mildly dependent on rriii also in these cases. In Table |^, we list the numerical values of 
the K-i&ctoTS for all the channels corresponding to the LHC ^/s = 14 TeV, 7 TeV and 
Tevatron ^/s = 1.96 TeV, obtained at rriii = 2mt for /i/j = ^p = I'^f^t with k = 0.5, 1, 2. 
The values of the -fC-factors for k = 1 will be used in the following. 

We check consistency of our iT-factor normalization in the large rriii region, by com- 
paring our prediction for the color-summed ti invariant-mass distribution with the NLO 
prediction. In Fig. ^ we plot the ratio of these two cross sections. In the former distribu- 

(c) 

tion, we include the if-factors, while we do not include the non-resonant diagrams A^V 
[c.f. (ii) of Fig. 0]. The latter distribution is computed for the on-shell ti productions. 



by MCONLO Q, |5[ with CTEQ6M PDFs with a scale choice ^p = hr = Jmf +pf.j. 
As can be seen, both cross sections are mutually consistent within 2% accuracy up to 
ma = 800 GeV. 

Including non-resonant diagrams in a way that the gauge cancellation holds effectively. 



^In the gluon radiation functions, the terms enhanced by plus-distributions or delta-functions as z — ;■ 1 
are common in Mf and pil], while non-enhanced terms differ. 

"in pit], effects of resummation of threshold logs are also examined and found to enhance the normal- 
ization at 10% level. See also |||, ||, |5|. 
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LHC 14 TeV 


LHC 7 TeV 




Tevatron 




K 


gglc=l] 


ggic=Sl 


qq 


gg[c=l] 


gg[c=Sl 


qq 


^^[-i] 


ggic=Sl 


qq 


0.5 
1 
2 


0.79 
1.14 
1.48 


1.02 
1.39 
1.75 


0.88 
1.16 
1.42 


0.87 
1.31 
1.75 


1.13 
1.60 
2.07 


0.89 
1.18 
1.45 


1.30 
2.11 
2.95 


1.72 
2.60 
3.48 


0.87 
1.18 
1.45 



Table 1: X-factor normalization constant for each channel {gg color-singlet, octet and qq) for the 
LHC y/s = 14 TeV, 7 TeV and the Tevatron, with setting the factorization and renormalization 
scales to /i/j — iip ~ nrrit with k = 0.5, 1, 2. 
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Figure 9: Ratio of the two color-summed tt invariant-mass distributions: The former is our 
prediction including the iT-factors but omitting the non-resonant diagrams, and the latter is the 
NLO prediction by MCONLO for on-shell top quarks. 



our final formula for the matrix element at parton level reads 



M 



bW+bW- 



XI ^ hW+hW) = VK 



■' tt -^ 'nr 



(2.25) 



<ic) 



lie) 



Ml^' and M)^J are given in Eqs. (|2]23|) and (|]2|) 



3. Event Generation and Top-Quark Distributions 

In this section, we present numerical evaluations of various kinematical distributions of 
the top-quark computed from the pp — )■ bW~^bW~ cross section, using the theoretical 
framework explained in the previous section. In particular we study the bound-state effects 
on these distributions. 

Our numerical calculations are carried out based on the MadGraph output |5^ which 
makes use of the HELAS subroutines [p^] for helicity-amplitude calculations. The original 
MadGraph output code has been modified to implement the color-decomposition and to 
include the bound-state effects via the Green functions. For the convenience of the readers, 
we collect the formulas necessary for decomposing amplitudes into the color-singlet and 
octet components of the ti (or bb) system in Appendix j^. In particular, we discuss how to 
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implement the color decomposition into the MadGraph notation. The bound-state correction 
factor 






c.f. Eq. (2.8), is pre-tabulated to save time for computing the momentum-space Green 
functions* . 

We perform phase-space integrations using BASES/SPRING |6C], or alternatively by 



adapting our code to MadEvent [p^, p^ utilities (ver. 4.4.42), where both tools are able 
to generate unweighted events at the partonic final-state level. For each event, we assign 
the specific color-flow according to an ordinary manner, except the color-singlet channel, 
as explained in Appendix ^ The generated events can be subsequently provided e.g. to 
PYTHIA for simulations of parton-showering and hadronizations. 

In the main body of this paper, we do not consider the decay of Ws but consider only 
the observables constructed from the bW'^bW'' final state. The VF-boson decays can be 
incorporated at the PYTHIA stage, where however the polarization of 1^-bosons cannot be 
taken into account. Alternatively, one can calculate the helicity amplitudes including the 
decay of W^-bosons by specifying a decay mode for each T^-boson. In Appendix ^, as a 
sample case, we examine the distributions of dileptons in the dilepton mode, where both 
W^s decay leptonically, and study the effects of VF-boson polarization and bound-state 
corrections. 

Below we show the results at the partonic bW~^bW~ final-state level. We do not discuss 
the parton-showering and hadronization effects, in order to concentrate on the examina- 
tion of bound-state effects. For the parton distribution functions, we use the CTEQ6L1 
parameterization with the LO evolution (1-loop running) of the QCD coupling constant. 
We set the renormalization and factorization scales to /i/j = hf = rrit and incorporate the 



X-factors obtained in Sec. 2.3 to the cross sections in the individual channels. [The final 



formula for the matrix element is given by Eq. ( 2.251 ).] We set the top-quark pole-mass, the 



(tree-level) on-shell top-quark width and the strong coupling constant as mt = 173 GeV, 
Tt = 1.49 GeV and a^(M^) = 0.1298, respectively. 

In Fig. p!o|(a), we plot the ti invariant-mass distribution in pp — t- bW~^bW~ production 
at y/s = 14 TeV. The tt invariant-mass m^j is defined as the invariant-mass of the final 
bW^bW~ system. The green solid line represents the full result which includes the bound- 
state effects as well as the i^- factors, and the blue dashed line represents the Born- level 



result (the LO prediction in the conventional perturbative QCD approach). Fig. lC(b) 



shows a magnification of the same cross sections in the threshold region. As shown in 



|2C, 21|, theoretically the bound-state effects can be seen most clearly in the shape of the 
m^t distribution in the threshold region. One can see that the cross section is enhanced over 
the Born cross-section significantly by the bound-state effects, and there appears a broad 
peak below the threshold corresponding to the ^^o resonance state in the color-singlet tt 
channel. Far above the threshold, the bound-state effects disappear and the cross section 
approaches the Born-level distributions, up to the iT-factor normalization. 



*The S-wave Green function depends only on \p\ but not on the direction of the three-momenta. 
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Figure 10: ti invariant-mass distribution in pp -^ bW^bW~' at y/s = 14 TeV. Green solid line 
is our full prediction and blue dashed line is the Born-level prediction. The NLO ti production 
computed by MCSNLO is also plotted in red dots. Right figure is the magnification of the threshold 
region. 
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Figure 11: The same as Fig. |lOl but for the LHC ^/i == 7 TeV 



In the same figures, we also compare our prediction with the NLO rrijf distribution 
computed by MCSNLO H, H with CTEQ6M PDFs and the scale choice of fip = fiR = 



ra'. 



+ p'^^. The latter prediction includes the full NLO QCD corrections (but not the 
Coulomb resummation) for the on-shell ti productions; we switched on an option of MCONLO 
to incorporate off-shellness of the top-quarks effectively by re-weighting the cross section 



by skewed Breit-Wigner functions |61], so that the cross section is non-zero below the 
threshold. (However, non-resonant diagrams are not incorporated.) Below and near the 
threshold, our prediction is much larger than the MCONLO prediction, due to the bound-state 
formation. The two cross sections become approximately equal from around mjf ~ 370- 
380 GeV up to larger mn. Note that, in Fig. |^, the contributions from non-resonant 



diagrams are not included in our full prediction, whereas in Figs. 10 they are included. 
Integrating the distributions over m^j, the total cross section by our full (Born-level) cal- 
culation is estimated as a^^^r+iw- ~ ^^^ P^ (633 pb), while we obtain cr^f = 816 pb as the 
MCONLO prediction. 



The shape of the m^j distribution at the LHC 7 TeV, shown in Fig. 11, is similar to 
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Figure 12: Two-dimensional density plots of the bW^ and bW~ invariant-masses, for the events 
with mtt < 370 GeV at the LHC y^ — 14 TeV. Left figure (a) is the Born-level prediction and 
right figure (b) is our full result (including the bound-state corrections and X- factors). The mean 
value and the root-mean-square value displayed in each figure are calculated for the events within 
the frame of the figure. 

that for the LHC 14 TeV. The total cross sections are estimated to be 158 pb, 106 pb and 
146 pb by our full, Born- level calculations and MCONLO, respectively. 



Let us examine other distributions of the top quark. From Figs. |lO|and 11, it is obvious 
that the phase-space region, where the bound-state effects are important, corresponds to a 
rather limited portion of the full top-quark events produced at the LHC. Thus, in various 
distributions formed by the full events, the bound-state effects may well be negligible in 
practice. In order to examine the bound-state effects closely, in the following we consider 
the events restricted by m^f < 370 GeV (except where otherwise stated), instead of con- 
sidering the full events. They amount to about 9% (8%) of the full events according to 
our calculation with (without) the bound-state corrections at the LHC 14 TeV. Due to the 
large ti cross sections at the LHC, still a large number of such near-threshold events would 
be accumulated. In this paper, we do not discuss the important subject of how to measure 
m^i in real experiments, which requires detailed studies of errors and fake solutions; one 



may find them in earlier studies |15, 17 1. 

One observes a characteristic bound-state effect in the (bW~^)-(bW~) double-invariant- 



mass distribution. In Figs. 12, we show the density plots of the invariant-masses of the bW~^ 
and bW~ systems, given by (a) the Born-level prediction and (b) our full prediction. In 
each figure, the number of events is normalized to 100,000 in total, and the number of events 
per bin (0.2 GeVxO.2 GeV) is plotted with graded colors. The Born-level prediction (a) is 
essentially determined by the product of the Breit-Wigner functions, hence the distribution 
is almost reflection symmetric with respect to the on-shell lines (ph -|-p^+)^ = m^ and 



{Pb+Pw-f 



mr 



By contrast, the distribution by our full prediction (b) is not symmetric 



and biased towards the configuration, where one of t or i is on-shell and the other has an 
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mtt < 370 GeV 


Full events 


k 


both [%] 


either [%] 


both [%] 


either [%] 


1 


37.0 (41.1) 


87.0 (87.6) 


46.8 (47.7) 


90.3 (90.4) 


2 


55.5 (59.0) 


95.2 (95.2) 


67.6 (68.3) 


97.0 (97.0) 


3 


64.2 (66.3) 


97.3 (97.1) 


76.3 (76.8) 


98.5 (98.5) 


4 


69.2 (70.3) 


98.1 (98.0) 


81.0 (81.4) 


99.1 (99.1) 


5 


72.3 (72.9) 


98.6 (98.4) 


83.9 (84.2) 


99.4 (99.4) 



Table 2: A fraction of events which satisfy \mbw — rnt\ < kTt for both or cither of the bW 
invariant-masses. The events with tti^j < 370 GeV as well as the full events at the LHC 14 TeV are 
considered. In the bracket is shown the result in Born-level. 



invariant-mass smaller than mj. In fact, such a configuration is known to be the dominant 
configuration just below the threshold in e~^e~ 



tt p^, 24 1 , although in that case deviation 



from the double Breit-Wigner distribution [Fig. |l2|(a)] is more prominent. (Note that, below 
the threshold, t and t cannot become simultaneously on-shell.) 

In order to quantify the correlated deformation of the double-invariant-mass distribu- 
tion, we count the fraction of the events for which both or either of the bW^ and bW~ 
invariant-masses satisfy 



\mbw - "T-tl < kVt 



(3.2) 



where k = 1,2, ...,5. These fractions are tabulated in Table ^ for our full prediction 
and for the Born-level prediction. The bound-state effect reduces the fraction for which 
both invariant-masses are close to on-shell more than the fraction for which either of the 
invariant-masses is close to on-shell. In the former case, the change of the fraction by 
the bound-state effect amounts up to about 4%. For comparison, we also tabulate the 
same fractions for the full events; in this case, the variation of the fractions are small 
and at most 1%. In any case, a proper understanding of this effect would be important, 
since it potentially biases the mass cut and may affect, for instance, the top-quark mass 
measurement with high accuracy. 

Let us explain the mechanism how the bound-state effects alter the double-invariant- 
mass distribution. As shown in Appendix A, the leading part of the tt amplitude A^jj has 
a form 



M^o. 



/ 1 1 



^/p^ 



+ 



mr 



P't 



^Jp^ + ml 



(3.3) 



where Pt = Pb ~^ Pw+ ^^^ Pi ^ Pi ~^ Pw- ^'^^ defined in the tt cm. frame, and p 
denotes the magnitude of the top-quark three-momentum in this frame. The first factor on 
the right-hand side {{m^ii—H)' 



r{c) 



withiJ = 2^p^ + m1 + V^l.^{r) [c.f. Eq. (|2l5|) ], denotes 



the Green function of the tt system. We suppress the top-quark width for simplicity. In 
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Figure 13: Top-quark momentum {p = \pt\) distributions in the partonic center-of-mass frame 
for a fixed m^j = 370 GeV, (a) for the cofor-singlet channel, and (b) for the octet channel. The 
solid lines represent our full prediction, after omitting the non-resonant diagrams. The dotted lines 
represent corresponding distributions at the Born-level. In Fig. (a), the distribution by the non- 
relativistic formula (replacing E' ^> E in the Green function) is also plotted (the red dashed line). 



the case that tt is in the singlet channel (c = 1), the potential energy between t and t is 
negative, VQQj-)(r) < 0. Therefore, the denominator of the Green function become close 
to zero (hence, the Green function is most enhanced) if p is somewhat larger than the 

i.e., p > pos- On the other hand, the second 



on-shell momentum pos = \ h^ff/^ 



m: 



t ' 



factor {p^ — \Jp^ -|- mf)^^ + {p^ — y^fP + rn^)~^ is most enhanced when p = pos> since 
Pt+P^ = mtf. Thus, there is a competition between the two factors on the right-hand 
side of Eq. ( |3.3| ). As a consequence, the dominant configuration is the one in which neither 
of the two factors are maximal. In fact, in the dominant configuration one of t and i is 
on-shell and the other is off-shell: p| = mf, p^ < mf, or, Pf < mf, p^ = ml. The effect is 
opposite in the case that ti is in the octet channel (c = 8). Since the magnitude of the octet 
potential is much smaller than the singlet potential, ^qcd/^qcd ~ ~l/8' the bound-state 
effect turns out to be much larger in the singlet channel than in the octet channel. 

Displayed in Figs. |l3|(a) and (b) are the top-quark momentum [p) distributions of 
the events with m^f = 370 GeV {not with mn < 370 GeV), for the color-singlet and octet 
channels, respectively. To see essential features, only the tt diagrams are taken into account 
and the i^-factors are not included. In each figure, the black solid (blue dotted) line shows 
our full prediction (Born-level prediction). The peak momentum for each distribution is 
shown with a vertical line. The peak momenta of the Born-level distributions are (to a 
good approximation) the on-shell momentum, Ppeak ~ POS = 65.5 GeV. We see that the 
bound-state effects shift the peak momentum by about 0.7 GeV to a larger value for the 
color-singlet distribution, while the peak momentum of the color-octet channel is shifted 
only by 50 MeV to a smaller value. In the color-summed cross section, the peak momentum 
is shifted to a larger value. Consequently, one of the invariant-masses of the bW~^ and bW~ 
systems is reduced below mt- The integral of this effect over the region TTi^f < 370 GeV 
can be seen in Fig. |T2|(b). 
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One may suspect that the above shifts of the invariant-masses (or the shift of the 
peak momentum) may be an artifact of our specific method to interpolate the ti cross 
sections in the threshold region and in the higher rritt region. To check this, let us estimate 
the size of the shift of the peak momentum at m^j = 370 GeV and compare it with the 
above prediction. The distance a top-quark propagates before it decays is estimated as 
7CT = mii/[2mtTt) = 0.72 GeV~^ = 1/(1.4 GeV). This distance is considered to be within 
the range where the potential VQ(jj-,(r) can be estimated perturbatively, although the 1- 
loop potential tends to underestimate the bound-state effect. ''^ The shift of the average 
momentum may be estimated by [{m^f — Vq(^j-,(7ct)}^/4 — ?n,f] — j»os = 0.66 GeV. 
Hence, the effect seen in Fig. [l3|(a) seems to be physical. 

We note that the effect elucidated here is a kind of effect that can never be seen in 
perturbative QCD computations for the on-shell ti productions, such as those given in 



1 34, ^, ^. This is because, the effect originates from the exchange of Coulomb gluons 
between off-shell t and on-shell i (or vice versa). Our full prediction correctly incorporates 
the (gauge-independent) LO off-shellness of the top quark as dictated by the exchange of 
Coulomb gluons, which is crucial for predicting the deformations of the top-quark momen- 
tum distribution and the double-invariant-mass distribution of the bW^ and bW~ systems. 
Now we are in a position to understand the origin of the abnormally large enhancement 
of the cross section, which we observed in Sec. 2.1, in the case that we use the non- 
relativistic formula for the differential cross section at large rritf, see the red dot-dashed 
line in Fig. S. The non-relativistic formula corresponds to replacing the Hamiltonian H = 
2i/p^ + "m-t + ^qcd('") by ^ = 2?Tii +p^ /mt + VQ'^j-,(r) in the Green function in Eq. (^). 
Thus, the non-relativistic formula overestimates the kinetic energy of the ti system in the 
large m^i region, 2mt +fP'/mt > 2\JfP^^rm^. For this reason, the two factors on the right- 
hand side of Eq. (^]^) can be brought close to maximal simultaneously with a nearly on-shell 
momentum, p ~ pos, since all the denominators in this expression nearly vanish. Since 
the individual factors are made of pole-type functions, applying an inaccurate kinematical 
relation only in one of the denominators can lead to a substantial overestimate of the cross 
section. In Fig. [T3|(a) we also plot the top-quark-momentum distribution computed with 
the non-relativistic formula, Eq. (2.8) after the replacement E' ^ E for the singlet channel 



(red dashed line). As can be seen, the peak momentum approaches the on-shell momentum 
and the distribution is more enhanced around the peak, compared to our full prediction. 

Other top-quark distributions are less affected by the bound-state effects. In Fig. p^(a), 
we show the normalized distribution of the top-quark momentum p = \pt\ = \pb + Pw+\ ™ 



the laboratory frame. In Fig. 14(b), we show the normalized distribution of the invariant- 



mass of bW {bW^ or bW ). The Born- level and full predictions are shown by the green solid 



and blue dashed lines, respectively. These lines in Fig. 14(b) correspond to the projections 
of Figs. |l2|(a,b) to the 'mf,^r+ (or miy^-) axis. All the histograms in Figs. [lj(a,b) are 
normalized, such that their integrals take the same value. 

In e~^e~ collisions, the top-quark momentum distribution in the threshold region is 



'See e.g. [fed] for the recent status of the perturbative prediction for the QCD potentiaL 
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Figure 14: Normalized distributions of (a) the top-quark momentuni, p = \pt\ = \pb +Pw+\i ^-nd 
(b) the bW invariant-mass, both defined in the lab. frame and for the events with rriti < 370 GeV. 
Green solid lines represent our full predictions, while the blue dashed lines represent the Born-level 
predictions. 

known to be proportional to the absolute square of the momentum-space Green func- 
|], whose shape is strongly influenced by the bound-state effects. At hadron 



tion [|5 



colliders, the top-quark momentum is boosted along the beam direction, and also the par- 
tonic collision energy is not fixed. As a result, even if we limit the events to those with 
"T-it < 370 GeV, the distribution of the top-quark momentum (p, defined in the lab. frame) 
is not much affected by the bound-state effects at hadron colliders. 

The mij\y distribution is important for the determination of the top-quark mass, hence 
it should be understood well. The bound-state effects deform the Born-level m^vv distribu- 
tion towards the lower side. The mean values of nibw over the range \mi,w — Tnt\ < 5 GeV 
are estimated to be 172.7 GeV and 172.9 GeV, for the full and Born-level predictions, 
respectively. The change of the mean value is about —200 MeV, for the restricted events 
with rritt < 370 GeV. At the LHC 7 TeV, we obtain almost the same result as in the 
14 TeV case. At the Tevatron 1.96 TeV, where the qq color-octet channel dominates, the 
mean values of mhw are estimated as 172.96 GeV and 172.98 GeV, respectively. Thus, the 
variation of the mean value is rather small. Note that MCSNLO predicts a mi,]y distribution 
similar to the Born-level distribution, since it simply re- weights the on-shell ti cross-section 
by skewed Breit-Wigner functions. 

4. Summary 

In the first part of this paper (Sec. 2), we explain our theoretical framework for including 
the bound-state effects in the fully differential cross sections for the top-quark production 
and their subsequent decay processes at hadron colliders. 

We formulate a theoretical basis to compute the fully differential top-quark cross- 
sections, which are valid at leading-order both in the threshold and high-energy regions, 
and which smoothly interpolate between the two regions. The tree-level ti double-resonant 
amplitude for each process is multiplied by a correction factor, which is written in terms 
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of the (well-known) non-relativistic Green function, but using a modified energy. This 
prescription preserves the required unitarity relation between the total and differential 
cross sections, which would be seriously violated had we used the naive non-relativistic 
formula for the differential cross sections at higher energies. 

We also include into the cross sections two important subleading corrections induced 
by the large top-quark width, (i) In addition to the tt double-resonant diagrams, which 
receives the above bound-state corrections, we include the contributions of non-resonant 
diagrams, whose effects are comparatively larger at higher energies. (This is more or less 
trivial.) (ii) As long as we perform numerical integrations of the differential cross sections, 
a sizable phase-space-suppression effect in the threshold region is inevitable, due to the 
sizable off-shellness of top quarks. In order to effectively account for the gauge cancellation 
by the Coulomb enhancement, we compensate the phase-space-suppression effect (by hand). 

Finally we incorporate ISR and the i^- factors. ISR effects are incorporated by con- 
necting our differential cross sections to a parton-shower simulator such as PYTHIA or 
HERWIG. We determine the iT-factors for the cross sections in the individual channels by 
matching them to the corresponding NLO ti invariant-mass distributions in the threshold 
region. With an appropriate choice of scales and using TTijf-independent i^-factors, we have 
checked that the color-summed ti invariant-mass distribution agrees with the conventional 
NLO prediction by MCONLO reasonably well at high energies. 

In the latter part of the paper (Sec. 3), using the above fully differential cross sections, 
we compute numerically various kinematical distributions of the top quark, constructed 
from the momenta of the bW~^bW~ final state (at the parton level). Our computations 
are carried out by MC event-generation using MadGraph, after implementing the color 
decomposition and the bound-state corrections to the output codes. 

We confirmed that our prediction reproduces the known NLO predictions for the ti 
invariant-mass distribution in the threshold region (which include bound-state effects) at 
the LHC 7 TeV or 14 TeV; in particular it exhibits the IS" resonance peak below thresh- 
old. Furthermore, our prediction approaches smoothly to the conventional NLO prediction 
(without bound-state effects) at higher invariant-masses, from around 30 GeV above the 
threshold. 

We restrict the events to those with ttt-jj < 370 GeV (in the case 2mt = 346 GeV), 
corresponding to about 10% of the full events, and examine kinematical distributions other 
than the m^j distribution. In particular, a characteristic bound-state effect on the (bW^)- 
{bW~) double-invariant-mass distribution is observed. The distribution is deformed from 
the double Breit-Wigner shape, towards the configuration in which one of the ti pair is 
close to on-shell and the other has a smaller invariant-mass than m^. 

The effect can be understood as a consequence of a competition between the contribu- 
tions from the (dominant) color-singlet Green function and from the t and i propagators. 
If the top-quark width were tiny, the Breit-Wigner distribution would tend to a delta- 
function, and the top quarks would be forced to on-shell. Due to the large decay width, 
however, the binding effect (towards off-shell mass) and the Breit-Wigner constraint (to- 
wards on-shell mass) remain to be competitive up to a few tens GeV above the threshold. 
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This effect lowers the mean value of each bW invariant-mass by a few hundred MeV for 
the above restricted set of events. The correlated deformation of the double-invariant-mass 
distribution may affect the mass cut and eventually the top quark mass measurement at 
the LHC. This requires further careful investigations. It would be worth emphasizing that 
the bound-state effect elucidated here can never be seen in the conventional perturbative 
QCD corrections to the on-shell ti productions, since the off-shellness of the top quark 
by the LO Coulomb binding effects plays a crucial role, and therefore it signifies a unique 
aspect of the present study. 

We examine other distributions, namely the (single) bW invariant-mass distribution 
and top-quark-momentum distribution. The bound-state effects on these distributions as 
a whole are not very significant, although there are certain systematic tendencies in the 
small deformations of the distributions, such as the aforementioned shift of the mean value 
of the invariant-mass. Furthermore, the dilepton distributions are examined including the 
leptonic decays of VF's at the matrix-element level (Appendix D). We have confirmed the 
previous observations that the effects of VF-boson polarization are quite significant, whereas 
we find that the bound-state effects are much smaller. 
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A. Green function of the ti system 

In this Appendix we explain how part of the Feynman amplitude corresponding to / — )■ 
ti — ;■ bW~^bW~ (/ = gg and qq) can be identified with a Green function that dictates 
the time evolution of the ti system. The argument is not restricted to the non-relativistic 
region. In order to define the bound-state as an eigenstate of the full Hamiltonian, we 
consider the limit F^ — t- 0.* 

Before providing a general argument, it would be pedagogical to demonstrate a decom- 
position of the three-point function (0 | T:ip{x)j^tp{x) : V'(y) V'(-2) | 0) in the free top-quark 
case, i.e., in the limit a^ — )■ 0. Here, ^p{x) represents the top-quark field operator. The 
three-point function is included, for example, in the amplitude for 7* — )• tt — )• bW^bW~ . 
One can express the three-point function as 

(0 \T.m-,-m:m «3) I 0),„ S f^-t"'- ," -f-'^:"'! (A.l) 

pf — mf + le pf — mf + le 



A+{p} A^ip) 



p^ — u}{p) + ie Pt + ^{p) — i^ 



7^ 



A_(p) A^ip) 



p^ — uj{p) + ie Pf + ^{p) 



le 



* We are not aware how to incorporate the resonance width in the formulation explained in this Appendix. 
It should not be a problem, however, since we want to find an expression that is valid as Tt -^ 0. 
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where F.T. stands for the Fourier transform, and io{p) = \/p'^ + rri^. In the second equahty, 
we have separated the contributions of particle propagating forward in time and antiparticle 
propagating backward in time. Here, 






represent the projection operators of a four-component spinor to the t and i components. 
Let us drop the contributions of non-resonant parts (far off-sheh contributions) in the last 
line of Eq. ( [A.l| ), which are 0{Tt/mt) had we retained the top-quark width in the top quark 
propagator. This corresponds to taking the contribution of the time ordering x^ < y^,z^ 
of the left-hand-side of the equation. Noting that p't +p\ = fnth the resonant part can be 
expressed as 

A.{p) _^^ K+{p) 



■r 



p^ — uip) A-ie ' Pf — uj{p) + ie 



mti-2uj{p) + i€ 



X 



+ 



p^ — uj{p) + ie Pf — oj{p) + ie 



(A.3) 



In coordinate space, this equation corresponds to the splitting of the time ordering x^ < 
iP , z^ into two orderings x^ < y^ < z^ and x^ < z^ < y^. On the right-hand side, the factor 
outside the square bracket represents the time evolution of the tt system. We may identify 
the denominator of this factor with m^j — Hq + ie and observe that the Hamiltonian Hq of 
the free ti system is given by 2uj{p) = 1\lp'^ + m^. The first term in the square bracket 
represents the propagation of t after t decayed first, whereas the second term represents 
the propagation of t after t decayed first. 

Now we develop a general argument. The four-point function of t and t 

(0 \T^{x{)^{x2)^{xz)i\){x^) I 0) (A.4) 

is a building block of the Feynman amplitude for I ^ tt ^ bW~^bW~ . The above four- 
point function can be decomposed to the sum of the different time orderings of xi, X2, X3, X4. 
As shown in |£^, the bound-state contributions are included in the orderings in which 
Min(x]',a;2) > yLax.{x^,x^). In fact, one may use the relation 



«(-°) = ^. (A-5) 
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to show that 

{0\Ti^{xi)ij{x2)i;{x3)i>{x4)\0] 



F.T. 



E 



1 



Mm{x1,x'^)>Ma.x{x'^,x1) 
^n{pf)(pn{pi) 



(A.6) 



2Mn mtt -Mn + ie 



1 



„0 

Ptj 



uj{pf) + ie 
1 



+ 



1 



P 



tj 



u}{pf) + ie 
1 



^L - w(Pi) 



le 



Ph 



uj{pi) + ie 



+ (non-resonant part), 



in the cm. frame of the tt pair, p^- and p^r {pj. and Pjf) denote the four momenta of 
the initial and final t (t), respectively, in the ti cm. frame; pi = pt^i = —pi^i and pf = 
Ptj = —Ptj- The first line on the right-hand side of Eq. ([A. 61) is identified with the Green 



function of the tt system, which includes the bound-state poles, {pf \ [mii — H + ie] 
The bound-state wave functions are defined as^^ 

(0|V(0)|t;p)(t;p|Vi(0)|n) 



'^n{p) 
^n{p) 



Mp) 

(n|V^(0)|t;p)(t;p|^(0)|0; 
2u:{p) 



(0 


V'(O) 


t] 


-p) {t; - 


-p\ 


^(0) 


n) 


(n 


V^(0) 


t; 


Mp) 
-p) {i; - 


-p 


V'(o) 


0) 



^ ' "• 2:(a '•- ■ ^^' 



\ t ;p) denotes the (anti)top-quark one-particle state with momentum p. The bound-state 
I n ) is defined as an eigenstate of the full Hamiltonian, H \n) = Mn | n ), and it is assumed 
to be a CP eigenstate. 

In Eq. ( |A.6| ), the bound-state poles stem from the time evolution between Min(x5,X2) 
and Max(x3,X4), whereas the single particle poles stem from the time evolution between 
Xi and X2 and between x^ and x^. For instance, in the case that x^ > X2 > x'^ > 
X4, one may insert the projection operators to the subspaces spanned by single particle 

s*^*^S' / i27vy'L{p) I *!^) ^^'P\ ^^^ I {iJ^Lip) I *! ~P) ^^5 -pI to extract the contributions 
from single-particle poles. Then the bound-state poles and single-particle poles appear 
from the Fourier transform 



{xl 



F.T. 



(t;p|V5(0) 



m* 



H + ie 



*^(-°--§)VS(o,f3)|t;-p) 
V5(o) I i- -p) 



(A.9) 



E 



1 (i;p| V(0)|n)(n|^(0)|t;-p; 



2M„ 



rrii 



Mn + ie 



+ (non-resonant part). 



^ifnip) is related to the Bethe-Salpeter wave function Xi(p) by 
'Pt -Pt 



X-n 



2 



d {x — y) exp 



,Pt-Pt, X , ,, ^° + y° 



(0|rV(x)i^(y)|n) 



p" — Lj{p) + ie p° — ijj{p) + ie 



+ (terms without a single particle pole). (A. 7) 
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and 

d{xl - x9,) (0 I V'(O) e-^^(^?-^°) I t;p) (A.IO) 

= -0 ^^r— -(0|V'(0)|t;p> 

+ (terms without a single particle pole). 

We may also express 9{x^ — x^) {i; —p\ e~*^'^3~^4) ^(0) | t;p) in a similar manner. The 
second line of Eq. ([A.9|) can be identified with the Green function of the tt system. 



B. Derivation of the off-shell suppression effect 

In this appendix we derive the off-shell suppression effect for the process I ^ tt with I = gg 
or I = qq in the threshold region. A similar formula for the process e^e~ — t- ti was derived 



in m 



The tree-level double-resonant amplitude has a form 

Mlf{pi,P2;Pb,Pw+^PhPW-)tree (B.l) 

= V{pt;pb,Pw+) ■SFiPt) ■'Pf\pi,p2]Pt,Pt) -SFiPt) ■'^(PuPb^Pw-)^ 

where V and V represent the decay part of t and i, respectively, Sp and Sp denote the 
propagators for t and i, respectively, in FWS: 



h + rnt) 
p'f -mf + imtVt 



Spipt) = -2 2 , ,^. -n ' (B.2a) 



pf -mf + imtVt 

(c) - 

Vj represents the tt production part which depends on the initial-state partons /. 

Integrating |A^Jr{ P over the bW~^bW~ phase-space for a fixed ii invariant-mass ^/s, 
we obtain the ti invariant-mass distribution in the following form: 

•^SfeHz+biy-^^) = J '^^P^'^'^P^-'^^ -^Sffls; ^t, Si), (B.3) 



where 



p{s) = 2^sT,{s)\Ap{s)\' , Af{s) = , (B.4) 

s — rrif + ZTTT-ti t 



and Tt{s) denotes the running width of the top quark with pf = s in the unitary gauge. 

2 

t 



Explicitly it is given by replacing mf by s in the on-shell decay width: 



3 



r,(.) = e{s-imw + m,f) ^A(l, ^, ^)/(^, ^) , (B.5) 



8\/2- 



vr 



where f{x,y) = l + x — 2y — 2x^ + xy + y'^ and A(a, b, c) = y^a^ + 6^ + c^ — 2{ab + bc + ca). 
In Eq. ( [B.3D , ^-^ og('^' **' •**) represents the off-shell tt cross-section, corresponding to the t 
and t invariant-masses of st and Sf, respectively. 
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The formula ( [B.3| ) is obtained in the following manner. First we define the t and t 
decay parts of the matrix-element squared, integrated over the bW phase-space, as 

^2\ 



J yPt 

Hk)^ /dcl>2(p?;P5,PH^-) J^(P^P) =^A-(l-75), (B.6b) 



where the summation is over the spins of the final 5-quark and VF-boson. Here and 
hereafter, ^4'' = ^^A'^'^^ denotes the Dirac conjugate. Decomposing the four-body phase- 
space and utilizing Eqs. (B.6), the bW~^bW~ cross-section is given by 



-£.M..EH.-(^) = il ^$^^^(«;^*'^*-)^ [{SUtSF} ■ Vf . {S,-a,Sl} . V\ 



(B.7) 



where Tr includes averaging over the spins and colors of the initial-state partons. The 
spinor matrices in the curly brackets {• • • } are calculated to be 

S\atSF = 2^Tt{pl)\AFipl)\^^St), (B.8a) 

Spa.S'j, = 27^r,(p?)|A^(p?)|2Sf(A-), (B.8b) 



^^.*-(^) = ^^^±-*±^^^- (^- 



with 



Thus, we obtain Eq. (|B.3| ) with the off-shell cross section 

^'i%is;plpf) = ^y"d$2(^;Pt,Pf)iv [s,(A) •T'f -EKA-) -^f '] • (B.io) 

Note that in the case Pt = P^ = fn^ "i off equals the on-shell ti production cross section. 

Below the ti threshold, the dominant kinematical configuration is such that either one 
of t and t is on-shell, because of the presence of the jAjT'p factors. Taking this into account, 
the off-shell cross section can be approximated by 



— —^ 'mi 771? 



<^itff («; su s,) ^ afUs) ^iil^l^ Fj(^, ^\ , (B.ll) 

^ Pt ^ St St ) 

where (Tj_^jj(s) are the cross sections for the on-shell ti productions; Fgg{x,y) = xy and 
Fgg{x, y) = {2 + X + y)/4. The factor A(l, ^, ^) originates from the ti phase-space volume 

and reduces to /3t = y 1 ^ in the on-shell limit st = St = ml . 

Thus, the ratio of the bW'^bW'^ cross section to the on-shell ti cross section is given 

by 



a 



(s) [dstds-t . . , ,A(l,f,f) ^(m\ m\ 



i^m+hw- - ' _ , -. -t p^st)p{st) --111^^-^ FjC-^,^]. (B.12) 



>) (5) J 27r2^^^ I3t '\st' Si 
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^/st^tist)^/si^i,{s^) in the p factors, which stems from the bW phase-space volumes, pos- 
sesses strong dependences on the invariant-masses of t and i: ^/sTt{s) ~ s^ since the 
running width behaves as Tt{s) ~ s^'"^. Below the threshold, either of t and t is forced to 
be off-shell (p^ < m^), hence y/stl't{st)^/sfTt{si) acts as a suppression factor. On the other 
hand, the remaining factors X//3t and Fj give only weak enhancement near the threshold. 
Thus, the effect of phase-space-suppression is approximately accounted for by the ratio 



VsiFtist) ^/s^^t{St) 



(B.13) 



rritTt mtVt 

for fixed Sj and sj. The ratio Eq. ( B.12| ) is considerably less than unity below the threshold. 



C. Color decomposition of the amplitude and color flow in gg channel 

In this section, we describe the color structure of the gg — >• ti matrix-element, and its 
decomposition into color-singlet and octet components. Moreover, we comment on the 
color flow in gg channel. 

By explicitly stating the color indices of initial-state gluons (a, b) and final-state t and 
i {i,j), the matrix element at the Born-level can be written in a following form; 



M^^iPk, Xk) = \ {t\ T']^^ Msipk, Xk) + \ [t\T' 



u 



MA{pk,Xk)-, 



(C.l) 



where Mg and Ma represent the subamplitudes for color-symmetric and anti-symmetric 
part, respectively, which depend on the 4- momenta and helicities of initial- and final-state 
particles. The color- factor in the color-symmetric part is decomposed into 



_ \ rpa rpb\ x<^b X _, _rabcrpc 



(C.2) 



where the first term represents color-singlet contribution and the second term color-octet. 
The color-anti-symmetric part contains only color-octet contribution; 



rpCl rpb 



_ pabcrpc 

2-' '^' 



(C.3) 



Since each part do not interfere with each other, the absolute square of the amplitude with 
summing over colors is given as a sum of each contribution; 



>J K 


2 

= E 


'}.N '^ 


colors ' ' 




+ E 


1 2 
jabcrpc 

2 '^ 



\MsY 



\ms? + y: 



_ rabcrpc 



\MaY 



(C.4) 



The first line in the r.h.s. in Eq. ( C.4 ) is the color-singlet contribution and the second line 
the color-octet. 
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^c; 2 




Figure 15: Color-flow diagrams in the gg — >■ ti amplitudes. Left diagram representing the color- 
factor S"'^6ij is for the color-singlet case, middle and right diagrams representing the color-factor 
{T°-T'')ij and {T''T"-)ij, respectively, are for the color-octet case. 



Alternatively, we may express the amplitude in the following basis; 



(C.5) 



The two basis, Eq. (^) and Eq. (^) are related by Ms = Mj^ + Mj^ and Ma 
Mj^ — Mj^. The absolute square of the amplitude is then written in a matrix form as,"f 



EKr=H 



colors 



\[MXMj, 




(C.6) 



By rewriting the absolute square of the color-singlet part of the amplitude in this basis, 
the corresponding matrix for the color-singlet gg — )■ ti amplitude is found to be 



E 



1 



S^^di 



\Ms[ 



3 ; ^x M-j^ 




(C.7) 



The color-octet contribution is obtained by subtracting Eq. (|C.7|) from Eq. (|C.6|) 



Finally, we comment on the color flow in gg channel. § In the color-singlet case, the 
color flow is disconnected between initial-state and final-state, reflecting the color-factor 



6 dij] see the left diagram in Fig. 15 



In the color-octet case, there exist two kind of color flows, middle and right diagrams 
in Fig. ^, associated with the color- factor; {T"-T^)ij and {T^T°-)ij, respectively. Either of 
the two may be selected according to the ratio |-/Vfj^p : l-Mjjp at each phase-space point 
and given helicities of initial and final-state particles in event generations. The color flow 
in qq channel is unique at Born-level. 

D. Dilepton distributions in dilepton decay mode 

In the main body of this paper we assume that the W^-bosons from top-quark decays are on- 
shell. In this appendix, we take into account decays of the W^-bosons at the matrix-element 
level. The advantage is to correctly take into account the off-shellness and polarization of 
the T^-bosons. This is crucial to predict correctly the angular distributions as well as other 



'"This matrix corresponds to the matrix CF in the MadGraph code (matrix. f). 

^We thank F. Maltoni for pointing out the modification of the color-flow structure in MadEvent, see also 
Ref. ||. 
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Figure 16: Distributions in an invariant-mass (left), a distance in rj-cj) plane (middle), and a 
difference in the azimutlial angle (right) of the two leptons for the dilepton events at the LHC 
y/s = 14 TeV. Top three graphs are for the full events, and bottom three graphs are for the events 
with mti < 370 GeV. Blue solid line is the full result, blue dashed line is the Born- level result, and 
red dot-dashed line is obtained by the full calculation but assuming unpolarized H^-bosons. 



kinematical distributions of the decay daughters of VF's. As an example, we show some 
distributions in the dilepton decay mode pp — )• bW~^bW~ — )• hi^v^hi'vi, where i = e, fi. 

We calculate the amplitudes from the set of Feynman diagrams obtained by just adding 
the W — )• ivi vertex to the Feynman diagrams for the bW~^bW~ final state. We generate 
the dilepton events at the LHC y/s = 14 TeV with standard kinematical cuts for the 
lepton momenta, \rji\ < 2.5 and px/ > 10 GeV. To avoid a singularity due to the vanishing 
running top-quark width, see Eq. ([B.5| ) in Appendix^, we restrict the phase-space integral 
region to rrihw > nib + mw ■ Note that the matrix-elements for rubw < rrib + mw are very 
suppressed. We set the VF-boson decay-width to Tw = 2.05 GeV. 



In Figs. 16, we plot three different distributions of kinematical variables constructed 
from the four-momenta of dileptons: (al, a2) the invariant-mass of the two leptons, (bl, 
b2) the distance in the r]-(j) plane, AR^ = a /Ar/I^ -|- A(/>|^, and (cl, c2) the difference of 
the azimuthal angles, Acpu. The first three graphs (al, bl, cl) correspond to the events 
from all the m^^ region, while the last three (a2, b2, c2) to the events with m^f < 370 GeV. 
The solid lines represent our full prediction, and the dashed lines represent the Born- 
level predictions. To see the effects of the non-zero VF-boson polarization, we also plot 
the distributions computed from the bW^bW~ events followed by the leptonic decays of 
on-shell unpolarized H^-bosons. 

We find that, due to the bound-state corrections, all three distributions are shifted 
to the lower side, although the variations are fairly small even for the events with m^f < 
370 GeV. By contrast, the effects of the non-zero VF-boson polarization is pronounced, 
especially for the events with rriu < 370 GeV. The most evident difference can be seen 
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in the lS.(j)u distribution, where the fuh calculation predicts that the number of events at 
^(t>u = is more than twice than the number of events at lS.(j)u = tTj while almost flat 
distribution by assuming the unpolarized T^-bosons. 



This finding is consistent with the similar study in Ref. [17| where the importance 
of the top-quark spin correlation in the ^(pu distribution is examined for the events with 
rritt ^ 400 GeV. They compare the A.(f)u distribution fully taking into account the top-quark 
spin correlation with that assuming the spherical top-quark decay into bW^s followed by the 
correlated VF-boson decay. Our calculation assuming unpolarized VK-boson decays includes 
the correct angular distributions in t — )■ hW decays, but forcing spherical distributions in 
W —7- Ivi decays. Thus, to predict the dilepton observables, all the spin correlations in 
decays of top-quarks and also VF-bosons are required. We confirm their finding that the 
difference in I^Ru distribution comes from mainly the difference in ^(pu distribution. 
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